Increasing Quantum Degeneracy by Heating a Superfluid 
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We consider a uniform superfluid confined in two compartments connected by a superleak and 
initially held at equal temperatures. If one of the two compartments is heated, a fraction of the 
superfluid will flow through the superleak. We show that, under certain thermodynamic conditions, 
the atoms flow from the hotter to the colder compartment, contrary to what happens in the fountain 
effect observed in superfluid Helium. This flow causes quantum degeneracy to increase in the colder 
compartment. In superfluid Helium, this novel thermomechanical effect takes place in the phonon 
regime of very low temperatures. In dilute quantum gases, it occurs at all temperatures below 
T c . The increase in quantum degeneracy reachable through the adiabatic displacement of the wall 
separating the two compartments is also discussed. 
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The thermomechanical effect is an important manifes- 
tation of superfluidity. It has historically been observed 
via the fountain effect, i.e. the increase in the pressure 
in a narrow tube, one of whose ends dips in a bath of 
superfluid Helium 4, when the tube is heated [if. It has 
recently also attracted interest in the context of dilute 
ultracold atomic gases @, as a potential signature of 
superfluidity in these systems. These works have focused 
on the dynamical aspects of the phenomenon, using the 
hydrodynamic Q or the classical field Q approach. 

The purpose of the present work is to investigate 
the thermomechanical effect by exploiting the conditions 
imposed by equilibrium thermodynamics, pointing out 
novel features exhibited by superfluids in properly cho- 
sen thermodynamic regimes. The experiment we pro- 
pose is reminiscent of the original fountain effect [ij, 
with one important difference. In the original experi- 
ment, the pressure at the surface of the liquid Helium 
bath is constantly equal to the saturated vapor pressure, 
and the height growth of the liquid in the narrow tube is 
determined by the equilibrium between gravity and the 
pressure increase Sp = sST of the liquid near the su- 
perleak. Here ST is the temperature difference between 
the bath and the tube, and s is the entropy per unit 
volume in the liquid phase. In the situation considered 
in this paper, the quantum fluid instead occupies a fixed 
total volume, and the flow of particles through the super- 
leak is caused by the density difference produced by the 
heating process. Consequently, it is determined by the 
compressibility of the fluid. In the case of liquid Helium, 
the compressibility is small; it is much larger in dilute 
quantum gases. We predict that, in the phonon regime 
of superfluid Helium, and for all temperatures below T c 
in the case of dilute gases, atoms flow through the su- 
perleak from the hotter to the colder region, contrary to 
what happens in the fountain effect, and resulting in an 
increase of quantum degeneracy in the colder compart- 
ment. The cooling mechanism proposed in the present 
paper, based on a filtering process through the superleak, 




FIG. 1. Schematics of the proposed experiment. The left and 
right compartments, initially in equilibrium, have constant 
volumes and are filled with a homogeneous quantum fluid. 
Heating the right compartment, or displacing the superleak, 
causes the system to evolve towards a new thermodynamic 
state satisfying chemical, but not thermal, equilibrium. 



differs from other adiabatic cooling mechanisms consid- 
ered in ultracold atomic gases, like the adiabatic forma- 
tion of Bose-Einstein condensation with non-harmonic 
traps ijj or the entropy exchange in mixtures of different 
atomic species [f|. 

We consider two compartments, hereafter called left 
(L) and right (i?) compartments, filled with a homo- 
geneous superfluid (liquid Helium 4 or a dilute atomic 
gas) and connected via a superleak, which only allows 
for the transmission of the superfluid component (see 
Fig. [H). Initially, the superfluids occupying the two com- 
partments have the same temperature (T° = Tjj = T°) 
and the same chemical potential (/i° = A*?j)- If the right 
compartment is heated, the system will eventually reach 
a new equilibrium configuration characterized by equal 
chemical potentials, but different temperatures. The 
equilibrium between the final-state chemical potentials is 
ensured by the flux of the superfluid through the super- 
leak. However, the condition of equal temperatures can- 
not be ensured because the superfluid component does 
not carry any entropy. 

By calling ST = Tr — Tl,Su = tir — tll and Sp = Pr — 
Pl the small differences between the final temperatures, 
densities and pressures of the two compartments, and 
imposing equal chemical potentials <5/i = fi(TR,riR) — 
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both partial derivatives being evaluated at the initial 
equilibrium parameters. The Gibbs-Duhem relation for 
a homogeneous system implies dp/dT\^ = s, where 
s = S/V > is the entropy per unit volume. As a con- 
sequence the pressure pn in the warmer compartment is 
always higher than pl . The density difference 5n can in- 
stead exhibit a different behavior. In fact, one can write: 
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the isothermal compressibility kt being a positive quan- 
tity. Hence, the direction of the flow is dictated by the 
sign of dp/dT\ n . If this derivative is negative, and thus 
dn/dT\ p > 0, the atoms flow from the colder to the 
warmer compartment (positive flow). The opposite hap- 
pens if dp/dT\ n > 0, and hence dn/dT\^ < (negative 
flow). In the latter case, the quantum degeneracy of the 
colder compartment increases. In fact, the total entropy 
in the left compartment remains constant, but the en- 
tropy per particle Sl/Nl decreases. In order to ensure 
adiabaticity, the heating process (as well as the displace- 
ment of the separating wall discussed in the last section) 
should be slow enough to ensure thermalization in the 
left compartment, and the velocity v of the atoms go- 
ing through the superleak during the process should be 
smaller than a critical value dictated by the geometry 
of the experiment. This condition will also ensure the 
absence of dissipation due to the creation of vortices. If 
the transmission of the atoms relies on quantum tunnel- 
ing, the current nv should be smaller than the critical 
Josephson current. The actual implementation of these 
conditions will be discussed elsewhere. 

Using the Gibbs-Duhem relation dp = sdT + ndu, 
and introducing the thermal expansion coefficient na p — 
—dn/dT\ p , Eq. ((2]) can also be written in the form 
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This form is particularly useful in the case of superfluid 
Helium where the quantities kt, s and a p are experimen- 
tally available. It explicitly shows that the direction of 
the flow depends on the sign of (kts — a p ). 

A general result concerning the behavior of the ther- 
modynamic quantity dji/dT\ entering Eq. ^ can be 
inferred in the very low temperature regime, where the 
thermodynamic behavior of a superfluid is governed by 
the thermal excitation of phonons (phonon regime). In 
this regime, the chemical potential behaves as p(n,T) = 
Ha{n) + ^ p hon(«, T) with 

7T 2 {k B Tf dco 
30 tfet dn 
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while no(n) and Co(n) are the T = values of the 
chemical potential and the sound velocity, respectively. 
In all known superfluids the quantity dco/dn\T is pos- 
itive. Therefore, in the phonon regime one always has 
du/dT\ n > 0, and hence a negative flow. 

We now discuss in a more systematic way the temper- 
ature dependence of dn/dT\^ in the case of superfluid 
Helium and of dilute quantum gases. 

Superfluid Helium. When one increases the tem- 
perature and leaves the phonon regime, the thermody- 
namic behavior of superfluid helium is soon dominated 
by the thermal excitation of rotons which deeply affect 
the behavior of the chemical potential. In Fig. [5] (left) 
we show the temperature dependence of fJ>(T) extracted 
from the measured thermodynamic functions of Helium 
Q for temperatures below 1.5 K. Its variation from the 
T = value at zero pressure, Hq/Ub = — 7.16 K, is small. 
However, its change of behavior when T increases, and 
in particular the change of sign of 8p,/dT\ n , is appar- 
ent. This effect is caused by the roton contribution to 
thermodynamics which can be written in the form 
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In Eq. © we have defined / = 2p\^fm* '/ (2irh 2 ) 3 / 2 , where 
Po is the momentum corresponding to the roton minimum 
in the excitation spectrum and m* is the correspond- 
ing effective mass [7J. The coefficient / depends only 
weakly on n 0] and, for simplicity, in Eq. ([S]) we have 
only retained the density dependence of the roton gap 
A(n) which is well known to decrease with increasing n, 
thereby providing a negative contribution to dfi/dT\ n . 
The approximate expression (/x p hon + Mrot) correctly re- 
produces the qualitative behaviour of the full thermo- 
dynamic result (see Fig. [2] left), which confirms the key 
physical role played by rotons. The curves refer to the 
value p = mu e n = 145.3 kg m~ 3 for the density, mue be- 
ing the the mass of a single Helium atom. The analysis re- 
veals that, for temperatures larger than ~ 1 K, superfluid 
Helium is charaterized by dn/ dT\^ > 0, i.e. by a positive 
flow through the superleak. In this regime, the thermal 
expansion coefficient is actually negative Q and the two 
terms in Eq. ^ add up with the same sign. Using the 
values of the thermodynamical functions given in [9l[loj|. 
we predict that the relative amplitude (T / p)dp/ dT\ ll is 
of the order of 1(T 2 for T = 1.8K. 

For smaller temperatures, when we approach the 
phonon regime, Helium exhibits a negative flow throw the 
superleak. The temperature where this effect is largest 
depends on the density, and decreases as the density in- 
creases. For the density p = 145.3 kg/m 3 , close to the 
lowest density for which Helium 4 is liquid, the relative 
amplitude of the predicted effect is largest for T = 0.8 K 
(see Fig. BJ where we find {T/p)dp/dT\ fl w -2 ■ 10~ 4 . 
This value is two orders of magnitude smaller than the 
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Temperature T [K] 

FIG. 2. Chemical potential as a function of T for three different superfluids. Left: superfluid 4 He, calculated using the 
equation of state in [g] (red). The orange curve is the sum of the phonon and roton contributions (Eqs. Q and ([5}); the 
phonon contribution is also shown separately (green). Middle: ultracold Bose gas (red), calculated from the Hartee-Fock 
prediction (Eq. JBJ) with gnjksTa = 0.3. Right: ultracold Fermi gas at the unitary limit (red) 11]. For both gaseous 
superfluids, the ideal gas prediction is also shown (green). 



typical relative amplitude corresponding to thermome- 
chanical effect in the roton region. The reason for this 
large difference can easily be understood from Eq. ([3]): 
in the phonon regime, kts and — a p have opposite signs, 
leading to a strong suppression of the total amplitude. 
Furthermore, both s and a p become smaller and smaller 
as T — » 0. Although small, the novel thermomechanical 
effect in the phonon regime of Helium should be large 
enough to be detected experimentally by measuring the 
change in the index of refraction or the dielectric con- 
stant e of the liquid in the left compartment. Indeed, 
it should lead to a relative variation Se/e ~ 10 -6 12(, 



which is well above the typical experimental sensitivity 
13]. It is worth pointing out that the thermomechanical 
effect discussed above is based on the assumption that 
the total volume Vl + Vr occupied by the fluid is kept 
constant. This effect should not be confused with the 
fountain effect, which occurs in situations where there 
is no constraint on the total volume. This latter effect 
is driven by the pressure difference sST and, hence, is 
always characterized by a positive flow through the su- 
perleak. 

Dilute quantum gases. The thermodynamic behav- 
ior of quantum gases is very different from the one of 
superfluid Helium, due to the absence of the roton min- 
imum in the excitation spectrum. As soon as one leaves 
the phonon regime described by Eq. (J4j) , the relevant exci- 
tations are of single-particle nature. In the case of dilute 
Bose-Einstein condensates, a reliable description of the 
single-particle regime is provided by Hartree-Fock the- 
ory, which yields the following expression for the chemical 
potential in the superfluid regime below T c 



fx(n, T) = g (n + 2n T ) = g 



C(3/2) 



mksT 



3/2 



(6) 

and the expression fi = [Md + 2gn for T > T c , /iid be- 
ing the chemical potential of the ideal gas. In Eq. ([5]), 
no and Ht are the densities of the condensate and of 



the thermal components respectively, g = A-Kh 2 a/m is 
the coupling constant (a is the scattering length char- 
acterizing the atom-atom interactions and m the mass 
of a single atom), and C(3/2) ~ 2.612. Consistently 
with the weak coupling scheme, in the second equal- 
ity of Eq. © we have used the ideal gas expression 
n T = C(3/2)(TOfc B r/27r?i 2 ) 3 / 2 for the thermal density 
[TBI ] . Equation ([6]) shows that the chemical potential is 
an increasing function of temperature up to the criti- 
cal temperature mksTc — 27r7i 2 (n/C(3/2)) 2 / 3 , where a 
reaches its maximum which is twice the T = value [16j ] . 
The temperature dependence of the chemical potential 
at fixed density can now be measured in quantum gases, 
where the equation of state is obtained by suitable inte- 
gration of the density profiles of the trapped gas or direct 
in situ measurements. These measurements have so far 
been carried out for Fermi gases (see FigJ5 (right)) [ill ] 
but can also be performed in Bose gases [17 1. 

Combining Eqs. <j2j) and (J6)), we obtain the following 
expression for the relative amplitude of the thermome- 
chanical effect in the Bose gas below T c : 



T dn 
n df 
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yielding the negative flow of particles through the super- 
leak for all temperatures [18|. Equation ([7]) also shows 
that the relative amplitude is of order 1 for tempera- 
tures of the order of the critical temperature. The flow 
through the superleak being negative in the whole range 
of temperatures below T Cl and the relative amplitude of 
the effect being very large, are the two main differences 
exhibited by dilute BEC gases with respect to superfluid 
Helium. In particular, the second feature is the conse- 
quence of the large compressibility of the gas. 

It is interesting to discuss now the thermomechanical 
effect in terms of the increase of quantum degeneracy in 
the left compartment. Although our expression (|6|) for 
the chemical potential accounts for the weak interaction 
between the atoms in the condensate, the entropy of the 
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FIG. 3. Amplitude ^ | of the thermomechanical effect 

in liquid 4 He: dependence on the temperature T for the 
fixed density p = 145.3 kg/m 3 . 
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gas in each compartment is well described by that of an 
ideal gas. In this approximation, the entropy of a Bosc- 
condensed gas depends on the temperature and the vol- 
ume but not on the number of particles. Therefore, the 
flux of particles through the superleak, caused by heat- 
ing the right compartment, changes neither the entropy 
Sl nor the temperature Tl of the left compartment. 
However, as Sl is distributed among a greater number 
of particles, the quantum degeneracy of the left com- 
partment increases. Assuming that Vr 3> Vl, we find 
S(S L /N L )/(S L /N L ) = -3/2(T a /T?f/ 2 5T R /T a , where 
Tj? is the critical temperature for the initial density uq. 



We now discuss an even more efficient way to increase 
quantum degeneracy by exploiting the novel thermome- 
chanical effect described above. We consider the adia- 
batic displacement of the wall separating the two com- 
partments of Fig. Q] towards the right. This displacement 
has two important consequences. First, due to adiabatic- 
ity, the temperature decreases in the left compartment 
and increases in the right compartment, the entropy of 
the gas being proportional to VT 3 ^ 2 if the effect of inter- 
actions can be neglected. Second, the number of atoms 
in the left compartment increases, due to the presence 
of the superleak. Note that the temperature variation 
in each compartment is independent of the presence of 
the superleak as the entropy does not depend on the 
number of particles. In the absence of the superleak, 
n,L would decrease, and the quantum degeneracy in both 
compartments would remain unchanged. In the presence 
of the superleak, the situation is different: we find that 
the density in the left compartment actually increases 
with respect to the initial value (in contrast to the de- 
crease one would naively expect due to the increase of 
Vl). This follows from the thermal density being higher 
in the right compartment and the chemical equilibrium 
condition fiL = /ir. These two combined effects cause 
a strong increase of quantum degeneracy. We express 
it in terms of the reduced temperature tl = Tl/T c l in 
the left compartment, where the critical temperature T c i 
corresponds to the final value of ul- Assuming that Vr 
always remains much larger than Vl , the thermodynamic 
properties of the right compartment are not affected by 
the displacement of the superleak, and the final value of 



tl is given by: 
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where tq is the initial reduced temperature, and V® and 
Vl are the initial and final volumes of the left com- 
partment, respectively. Hence, the quantum degener- 
acy can in principle be made arbitrarily large by taking 
Vr » V L » Vl 

The novel thermomechanical effects we have predicted 
for the Bose gas are also expected to occur in a Fermi gas 
in the superfluid regime. The most interesting case is the 
unitary Fermi regime, corresponding to an infinite value 
of the scattering length. This system has been achieved 
experimentally tllj and its theoretical properties are rea- 
sonably well understood An interesting feature of 
this system is that its thermodynamic functions exhibit 
a universal behavior, independent of the value of the in- 
teraction strength. In Fig. [5] (right) we report the recent 
experimental measurement of the chemical potential as 
a a function of T/T F , where k B T F = {h 2 /2m){3ir 2 n) 2 / 3 
is the Fermi energy defined in uniform matter, together 
with the prediction for the ideal Fermi gas. Like in the 
case of Bose gases, the chemical potential of the Fermi su- 
perfluid is an increasing function of T up to the critical 
value. Its measured behavior is reasonably well repro- 
duced by the theoretical prediction obtained using the 
variational many-body formalism [l9| • In the case of the 
ideal gas, the chemical potential is instead a decreasing 
function of T, revealing that the positive slope of fi(T) 
is a clear consequence of superfluidity. Compared to the 
Bose case, the relative amplitude {T/n)dn/dT\^ of the 
thermomechanical effect is smaller in Fermi gases because 
of the much smaller compressibility of these systems. 

Finally, we briefly discuss how the thermomechanical 
effect could be implemented in a trapped quantum gas. 
Well established tools, such as magnetic traps and optical 
dipole potentials, can readily be used to create the two- 
compartment potential 1 20] , including the implementa- 
tion of box-like traps [2l[. The superleak can be im- 
plemented using quasi-ID optical dipole potentials with 
a reduced number of transverse excitation modes [22| . 
Other important issues concern the consequences of the 
thermomechanical effect in the presence of harmonic con- 
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finement, and in particular the optimization of the cool- 
ing mechanism through the adiabatic displacement of the 
wall separating the two compartments. 
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